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ABSTRACT 

We investigate the impact of nonaxisymmetric structure on estimates of galaxy inclinations 
and position angles. A new minimization technique is used to obtain estimates of inclination 
and position angle from a global fit to either photometric or kinematic data. We discuss possible 
systematic uncertainties which are much larger than statistical uncertainties. Our investiga¬ 
tion reveals that systematic uncertainties associated with fitting photometric data dominate the 
formal statistical uncertainties. For our sample of inclined galaxies, we estimate that nonaxisym¬ 
metric features introduce inclination and position angle uncertainties of « 5°, on average. The 
magnitudes of these uncertainties weaken the arguments for intrinsically elliptical galaxy disks. 


Subject headings: galaxiesdundamental parameters 
namics 

1. Introduction 

We observe disk galaxies projected at random 
angles on the sky and must, therefore, determine 
the inclination geometry for each galaxy in order 
to deproject line-of-sight velocity data, or discover 
the intrinsic shape of the disk, etc. The in-plane 
orbital speed is needed for mass estimation, and is 
also one of the parameters required for the Tully- 
Fisher relation (hereafter TFR, Tully & Fisher 
1977). Furthermore, correction for internal extinc¬ 
tion is also inclination-dependent (e.g., Giovanelli 
et al. 1994), affecting the total luminosity of the 
galaxy, the other parameter of the TFR. Thus er¬ 
rors in both the inclination and position angle, 
propagate directly into dynamical mass estimates, 
luminosity estimates, and contribute to the scatter 
in systematic relations. 

The most common procedure to determine ori¬ 
entation geometry is to fit ellipses to a photomet¬ 
ric image of the galaxy, and to infer the projec¬ 
tion geometry by assuming the disk to be intrinsi¬ 
cally thin and axisymmetric. The ellipticity of the 
projected image, e = 1 — b/a where b and a are 
respectively the semi-minor and semi-major axis 
lengths, is related to the inclination angle through 


galaxies:photometry — galaxies:kinematics and dy- 

cosi = 1 — e, for an assumed intrinsically axisym¬ 
metric, infinitesimally thin, disk. Here, i is the 
inclination angle of the disk plane to the plane 
of the sky {i = 0° for a face-on galaxy), and we 
use the symbol (p for position angle of the line of 
intersection of these two planes. A slightly more 
complicated relation exists for galaxies that are as¬ 
sumed to be thin oblate spheroids (Hubble 1926). 

If 2-D velocity fields are available, they can also 
be used to determine inclinations and position an¬ 
gles. The previously mentioned assumptions of 
intrinsic axisymmetry and thinness are again typ¬ 
ical. A drawback to this method is that it is poorly 
suited to determine inclinations where the rotation 
curve rises slowly (Van Moorsel & Wells 1985). On 
the other hand, 2-D velocity fields can clearly de¬ 
lineate position angles. 

Garcia-Gomez & Athanassoula (1991) sum¬ 
marize several other procedures that have been 
adopted to estimate e from photometry in the 
presence of the complicating non-axisymmetric 
features - spirals, bars, lop-sidedness, etc. - man¬ 
ifested by virtually all disk galaxies. 
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The IRAF^ utility ellipse implements the 
method outlined by Jedrzejewski (1987) and fits 
ellipses to a number of isophotes independently, 
although the task can be instructed to hold any 
parameter, such as the photometric center, fixed 
for all. The formal statistical uncertainties due 
to noise in the data (hereafter abbreviated to sta¬ 
tistical uncertainties) in the fitted parameters are 
not easily obtained (Busko 1996), but are gen¬ 
erally quite small. The user of this utility then 
has to make a judgement as to which isophote, or 
group of isophotes, offers the best estimate of the 
projection geometry. A common prejudice is that 
the outer disk is more likely to be intrinsically 
round, leading to a preference for the outermost, 
or the average of several outer ellipses. Another 
strategy is to identify the true shape with a radial 
range over which the fitted parameters do not vary 
much. 

Such a procedure clearly does not utilize all 
the information in the image. More sophisti¬ 
cated methods include: minimizing the m = 2 
component of the Fourier transform of an image 
(Grosbpl 1985), bulge-disk decomposition (Kent 
1985), maximizing the m = 0 component of depro- 
jected images (lye et al. 1982), and fitting a para¬ 
metric model for the disk (Palunas & Williams 
2000). We propose another global method which, 
like those just cited, assumes that the galaxy disk 
is infinitesimally thin and intrinsically axisymmet- 
ric. Our method returns e and </> values from a 
chi-squared minimization using a large part of the 
image. The details are summarized in the Ap¬ 
pendix. 

The statistical uncertainty of the methods de¬ 
scribed above is usually small, but it is clear that 
the possible systematic error could be large if the 
disk is intrinsically non-circular or if nonaxisym- 
metric features are present. In this paper, we at¬ 
tempt to quantify the possible systematic error 
from a detailed re-analysis of the photometric and 
2-D kinematic maps of a sample of 74 spiral galax¬ 
ies already reported by (Palunas & Williams 2000) 
(hereafter PW). 

After briefly introducing the data sets used in 
the study in §2, we give a general discussion (§3) 
of the statistical and systematic errors that can 

^IRAF is distributed by NOAO, which is operated by 
AURA, under contract to the NSF. 


arise when estimating projection angles. In §4, we 
present the results of photometric and kinematic 
fitting procedures applied to the PW galaxies. Un¬ 
like these authors, however, we study possible sys¬ 
tematic errors in estimating the projection param¬ 
eters by comparing our fits to the kinematic and 
photometric data. 

2. Data 

We have used the disk fitting methods described 
in the Appendix to estimate the projection angles, 
i and (j), from the /-band photometry and 2-D ve¬ 
locity fields for the 74 galaxies in the PW sample. 
PW kindly provided the reduced photometric im¬ 
ages with the foreground stars removed. 

PW also determined 2-D velocity fields for the 
same galaxies from Ha spectroscopy taken with 
the Rutgers Fabry-Perot imaging spectrophotome¬ 
ter. They obtained several (8 to 15) images spaced 
by 1 A(« 50 km/s at Ha) of each galaxy from 
which they derived data cubes (x,y,u;os)- These 
cubes were then reduced to 2D spatial maps by 
fitting the line profiles at each pixel with a Voigt 
function. Again, PW kindly supplied us with the 
resulting velocity fields, Ha intensities, continuum 
intensities, and their respective uncertainties. 

Inclined galaxies were preferred in this sample, 
which is therefore notably deficient in nearly face- 
on systems. This selection bias makes our study 
complementary to those by Rix & Zaritsky (1995) 
and Anderson et al. (2001) who focused on face-on 
systems. 

3. Statistical &: Systematic Uncertainties 

Noise in the data (photometric or kinematic) 
creates statistical parameter uncertainties, ax, 
which can be measured in a variety of ways. For 
example, the curvature of the reduced chi-squared, 
Xr , surface at the minimum, Xr = Xrmitu 
better, the distance from the minimum at which 
Xr = Xrmin + allows estimates of the sta¬ 
tistical uncertainies in the parameters, and any 
covariances. Confidence limits can be determined 
from the contours at = njv, where v is the 
number of degrees of freedom, i.e., the number of 
pixels minus the number of free parameters in the 
fit, and n is an integer. Clearly, the magnitudes 
of the statistical uncertainties in the estimated 
parameters are small when many pixels are used. 
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and statistical uncertainties in the projection an¬ 
gles frequently turn out to be ^ 1°. 

Sources of systematic uncertainties, Sx, in pho¬ 
tometric ellipse fitting include: (1) seeing, (2) 
dust obscuration, (3) finite disk thickness, (4) non- 
axisymmetric features in the disk, and (5) intrinsic 
ellipticity of the disk. Velocity field fitting is also 
subject to similar systematics. Points 4 and 5 are 
relevant to both methods. Additionally, if a nearly 
edge-on galaxy with finite thickness has rotational 
velocities that decrease with height away from the 
midplane, fitting a thin disk model leads to an ap¬ 
parent inclination that is somewhat less edge-on 
than in reality. 

Seeing tends to bias elliptical isophotes towards 
becoming round. While it is most severe in the 
central regions, the effect can be of the same order 
as an intrinsic ellipticity for moderate distances 
from the center (Trujillo et al. 2001). Using the 
typical seeing value quoted in PW of 1".5 and Fig¬ 
ure 1 of Trujillo et al. (2001), we find that at a ra¬ 
dius 12 times the seeing length a measured e value 
is « 90% of its actual value. This can lead to errors 
in f of a few degrees, depending on the value of i. 
However, since most of the pixels utilized to make 
the fits lie beyond 12 times the seeing length (and 
are less affected by seeing), the effect is probably 
less important than others we discuss below. 

Since inclined galaxies suffer more extinction 
than do face-on galaxies, dust is a source of sys¬ 
tematic error. Its effect can be reduced by utilizing 
photometry in near IR bands where dust extinc¬ 
tion is less severe, as we do in this study. 

Using simulation and analytical work, we have 
found that photometrically determined i values 
are significantly affected by finite thickness only 
for the most edge-on galaxies (inclination > 80°). 
Photometric fits of model inclined oblate spheroids 
agree well with the formula given in Hubble 
(1926) relating inclination, apparent axis ratio, 
and spheroidal axis ratio. Since the majority of 
the PW galaxies have i < 80°, we ignore this 
complication hereafter. 

Nonaxisymmetric structures in disks, such as 
bars, spirals, and warps, are a major source of 
systematic error. Stock (1955) investigated the 
impact of spiral structure on the photometric esti¬ 
mation of position angles, particularly the case of 
highly inclined galaxies, finding that spiral struc¬ 


ture can introduce uncertainties in (j) values of 
w 5°. We refer to the impact of all such nonax¬ 
isymmetric structure on the photometric projec¬ 
tion angles as the spiral effect. 

Several studies, most notably Binney & de 
Vaucouleurs (1981); Grosbpl (1985); Franx & de 
Zeeuw (1992); Rix & Zaritsky (1995); Beauvais & 
Bothun (1999); Anderson et al. (2001), have found 
hints that the disks of spirals are not perfectly ax- 
isymmetric, which may arise because the gravita¬ 
tional potential well within which the disk mate¬ 
rial orbits is intrinsically non-axisymmetric. We 
distinguish this global intrinsic oval distortion of 
the disk from more obvious features in the disk, 
such as spiral arms and bars. One reason for 
drawing this distinction is that the spiral effect 
can be investigated using photometric data only, 
while the effects of intrinsic ellipticity should be 
most clear when photometric data are compared 
to kinematic data (Franx & de Zeeuw 1992). De¬ 
pending on how the major axis of an elliptic disk is 
projected, an intrinsic ellipticity can lead to either 
larger or smaller inferred ellipticity values, leading 
to over- or under-estimates of the inclination, and 
to corresponding errors in the position angle. 

4. Results 

PW determined the projection angles for each 
galaxy separately from the photometric and kine¬ 
matic data. In general, they found somewhat dif¬ 
ferent values from the two methods, which was 
inconvenient for their objective. They wished to 
compare the observed rotation curve with that 
predicted from a photometric disk mass model, 
which requires a common projection geometry for 
both. They argued that the inclination and posi¬ 
tion of the geometric center could be determined 
most reliably from the photometry, while the po¬ 
sition angle of the major axis was best determined 
from the kinematics. They then re-fitted both 
data sets with the preferred parameter from the 
other set held fixed. 

Our inclinations derived from photometry are 
in good agreement with the values reported in 
PW. While reassuring, the agreement is scarcely 
surprising as we are using similar methods to fit 
the same photometric data. Discrepancies arise, in 
part, because we have fitted for both i and (j) while 
in PW <j) was held fixed at the value determined 
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from the kinematics, but the good agreement (the 
average difference between our Zphot values and 
those of PW is 2.9° shows that this constraint has 
a small effect. The largest difference is 19° and 
occurs for ESO 510G11, a galaxy with very strong 
spirals that extend into the outer parts of the disk. 

One of our principal objectives here is to un¬ 
derstand the discrepancy between the photomet¬ 
ric and kinematic values of the projection angles 
for each galaxy. The magnitudes of these differ¬ 
ences are illustrated in Figures la and lb. The 
statistical uncertainties are much smaller than the 
symbols in both figures. 

The inclinations show considerable spread 
around perfect agreement (Figure la). The aver¬ 
age photometric/kinematic inclination difference 
is « 5° with a standard deviation of similar mag¬ 
nitude. A comparison of position angles (Figure 
lb) shows a similar spread, with one extreme out¬ 
lier which we discuss briefly below. The average 
difference between position angle estimates is w 4° 
(excluding the extreme outlier). Note, however, 
that this galaxy sample is deficient in nearly face- 
on galaxies and that Anderson et al. (2001) found 
larger differences in a sample of 10 low inclination 
(z < 20°) galaxies. 

Figures 2a and 2b show the relations between 
the inclination and position angle differences and 
the photometric inclination of the galaxy. While 
there is significant spread in Figure 2b, there 
appears to be a trend for more highly inclined 
galaxies to have smaller position angle differences. 
Such a trend is expected, since any nonaxisym- 
metric structures present in highly-inclined galax¬ 
ies should have only a small effect on (/)phot- This 
trend also agrees with the aforementioned larger 
differences found for low inclination galaxies (An¬ 
derson et al. 2001). 

The largest position angle difference in Figure 
2b corresponds to a galaxy which has a low in¬ 
clination and strong spirals which extend to the 
outer reaches of the disk. It is a prime example of 
the kind of systematic error that can be attributed 
to the spiral effect. Additionally, this kind of mis¬ 
alignment could cause serious problems if single¬ 
slit spectroscopic data were to be utilized to deter¬ 
mine rotation curves. Misalignment between the 
slit and the true major axis results in an under¬ 
estimate of rotational speed and hence dynamical 
mass. 


4.1. Fixed Centers vs. Floating Centers 

In all the work reported elsewhere in this paper, 
we have required the rotation center of our fit¬ 
ted model to coincide with the peak of continuum 
brightness in the image (see Appendix). In this 
subsection, we report the consequences of relax¬ 
ing this requirement and determine the kinematic 
center from the velocity field alone. The offset be¬ 
tween the kinematic and photometric centers, and 
the resulting differences in estimates for z and (j) 
are illustrated in Figure 3. 

Figure 3a shows that the kinematic center is 
< 2."7 from the photometric center in the major¬ 
ity (67 of 74) of galaxies. The galaxy for which 
this distance is largest, ESO 322G42, has a slowly 
rising rotation curve, which creates nearly paral¬ 
lel isovelocity contours; the location of the cen¬ 
ter and the overall systemic velocity are poorly 
constrained by such data, but the inner slope 
of the rotation curve is independent of the cen¬ 
ter location. Six other galaxies are in the tail 
(Arcen > 2."7): One (ESO 445GI9) has a strongly 
asymmetric inner velocity field, while the others 
are either nearly edge-on and have roughly paral¬ 
lel isovels and slowly rising rotation curves, or lack 
Ha emission near the center; the location of the 
rotation center is not tightly constrained in the 
absence of kinematic information near the center. 

As evidenced by the histograms in Figures 3b 
and c, the projection angles that result from fitting 
with the kinematic center free to move agree with 
those from the fixed-center fits for the majority of 
the galaxies. In fact, none of the <j) values changes 
by more than 4°. The changes in inclination show 
a larger range, but 70 of 74 galaxies have Az < 7°. 
The galaxy with the largest Az shows evidence of 
being tidally distorted, having a long filament of 
anomalous velocities extending from its edge. The 
two galaxies with the next largest Az values both 
have strong spirals resulting in velocity fields with 
strong distortions. The remaining galaxy in the 
tail of the Az histogram is a highly inclined galaxy 
that is also in the tail of the Arcen histogram. 

We find that unless the kinematic center is 
poorly constrained by a slowly rising rotation 
curve, it lies within « 2."5 of the position of the 
photometric center for most galaxies. So, if rota¬ 
tion curves of our galaxy sample were determined 
from single-slit spectra taken through the point 
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of maximum continuum, they would not be sig¬ 
nificantly biased towards rising slowly. Allowing a 
small shift in the kinematic center generally makes 
little difference to the estimated projection angles, 
except in a few special cases. 

4.2. Spiral Effect Systematics 

In an attempt to estimate the magnitude of the 
uncertainties due to the spiral effect, we have pro¬ 
ceeded as follows. (Note that these values must 
be viewed as lower limits to the true systematics, 
which are unknown.) We use a first estimate of 
the photometric projection parameters from our 
Xr minimization procedure to deproject a galaxy 
image. We then rotate it through an angle (3, so 
that the nonaxisymmetric structure lies in a new 
orientation, and then reproject it to the original 
best-fit i and (j). We then apply our photometric 
fitting apparatus to this new image to find new 
best-fit i and (j) values. We interpret the spread in 
both i and </> as /3 varies through a half-rotation as 
the measures of systematic uncertainties Si and 
due to the presence of nonaxisymmetric features 
in the image. A few representative plots of i and 
(j) versus (3 are presented in the next section. 

The results of this investigation are shown in 
Figure 4. Galaxies that have i ^ 75° have been 
excluded from this analysis because deprojection 
must be considered suspect for such large inclina¬ 
tion angles. The general trend of decreasing Si and 

with increasing i concurs with our earlier asser¬ 
tion that nonaxisymmetric structure should com¬ 
promise the model fits more modestly for highly 
inclined galaxies. The galaxy with the highest val¬ 
ues of Si and (ESO 323G39) appears to have a 
significant bar. 

The numbers in Figure 4 denote different galaxy 
T-types, where low numbers indicate early-type 
spirals and high numbers indicate late-type spirals 
and irregulars. The different T-types do not nicely 
segregate in either frame of Figure 4. However, 
there is a statistical trend for later type galaxies 
to have larger Si and values. Early type spi¬ 
rals (T-types 0-3) have average Si values of « 3° 
and average values of « 4°. Spirals with T- 
types between 4 and 6 have average Si « 4° and 
« 4°. Late type spirals (T-types > 7) have 
average Si w 6 ° and s^ « 7°. The differences 
between the means of the early and intermediate 
spirals are not statistically significant. However, 


the differences between the early and late types’ 
means are significant at the 95% level. For all 
galaxies with i < 75°, the average Si « 4° and 

« 5°, with a large spread for both values. The 
magnitudes of Si and suggest that the uncer¬ 
tainties in i and (j) values derived from photometry 
are dominated by systematic uncertainty due to 
nonaxisymmetric structure. 

4.3. Intrinsic Ellipticities vs. the Spiral Ef¬ 
fect 

We find differences between the projection an¬ 
gles determined by fitting photometric and kine¬ 
matic data. Franx & de Zeeuw (1992), Anderson 
et al. ( 2001 ), and others suggest that such dis¬ 
agreements are evidence of intrinsically elliptical 
disks, but such a conclusion requires that the pro¬ 
jection angles can be determined with sufficient 
precision for their differences to have meaning. 

Our investigation of the spiral effect points to 
significant systematic errors in the photometric 
determination of projection angles which are of 
the same magnitude as the errors suggested by 
Franx & de Zeeuw (1992) to explain the scatter 
in the TFR. This leads us to wonder whether dif¬ 
ferences between photometric and kinematic pro¬ 
jection angles are evidence of anything more than 
our ignorance of the true angles. 

Fortunately, we can use our 2D velocity fields to 
search for evidence of intrinsic ellipticity. Franx, 
van Gorkom, & de Zeeuw (1994) formulate ellipti¬ 
cal orbits in disks in the epicyclic approximation 
with a small m = 2 forcing term added to the 
central attaction. This distortion is entirely con¬ 
tained within the plane of the disk. For the special 
case of a fiat rotation curve, the distortions in the 
radial and tangential directions are equal and can 
be described by, 

Vn = i^circCint sin 2 (V> - 7 ), (la) 

and 

= Wcirc[l + Cint COS 2 (V’ - 7)]- (lb) 

Here r^circ is the circular speed, Cint is the intrinsic 
ellipticity of the disk, tp is the angle in the disk 
measured from the projected major axis, and 7 is 
the angle between the projected major axis and 
the elliptical disk minor axis. When such a disk is 
inclined, the line-of-sight velocity field appears as 
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that of an axisymmetric disk with small correction 
terms added, 

'^los,ellip — 

^los,axisym 

Eint^^circ sini(cos 27 COS f/' — sin 27 sin ^), (2) 

where uios,axisym is the line-of-sight velocity for an 
axisymmetric disk with the same projection an¬ 
gles. Equation 2 demonstrates that a mildly el¬ 
liptical disk creates an m = 1 distortion in the 
residual velocity field, at a position angle which 
depends on 7 . (There is also an m = 3 term in 
the more general case of a non-flat rotation curve, 
see Franx, van Gorkom, & de Zeeuw 1994 for de¬ 
tails.) 

In order to quantify the strength of such dis¬ 
tortions, we introduce two quantities, P and Q 
which we derive from the velocity residuals in the 
following way. We divide the residual field inside 
the outermost fitting ellipse into 15 elliptical an¬ 
nuli, and divide each annulus into quadrants. We 
label the first quadrant (I) to be that bounded 
by the receding side of the projected major axis 
and the minor axis counter-clockwise from it, and 
the remaining quadrants (II, III, IV) are in the 
usual order. For each quadrant in each annulus, 
we calculate the average of the residuals fn^quad, 
normalized by the average of the magnitude of the 
residuals in the entire field. We then define P and 
Q values for each annulus as follows, 

P = \^n,l + rn,U — ^nJII + ^n.ivl (3a) 

and 

Q = \rn,l + TnJV - rn,ll + I'nJIll- (3b) 

Basically, P measures the strength of an m = 1 
distortion across the major axis, while Q measures 
the same thing across the minor axis. Plotting 
A = + values versus the semi-major axis 

of the annulus provides an asymmetry profile. 

Using equation 2, we have constructed a model 
residual velocity field for a disk with nested similar 
concentric elliptical orbits. The resulting asymme¬ 
try profile rises rapidly in the central region of the 
residual field, but then flattens and remains con¬ 
stant for the remaining radial range. None of the 
galaxies we have studied show this kind of asym¬ 
metry profile. Galaxies which have bars and/or 
strong spirals show peaks in their asymmetry pro¬ 
files (for example. Figures 5c VI and 5d VI), but 


the distortion is localized and not seen through¬ 
out the entire disk. It is also interesting to note 
that a higher-than-average level of asymmetry in 
a velocity field is not always associated with a sig¬ 
nificant difference between photometric and kine¬ 
matic projection angles. 

5. Individual Galaxies 

Here we present several individual galaxies that 
illustrate the difficulty of measuring accurate pro¬ 
jection angles. For each galaxy, we show an I- 
band image, a 2D velocity field, a residual pho¬ 
tometric image, a residual velocity field, a plot 
illustrating the spiral effect, and an asymmetry 
profile. The photometric images are presented in 
grayscales with units of magnitudes per arcsecond 
squared. The white circles in the photometric im¬ 
ages mark the locations of stars that have been 
removed. In the photometric residual images, the 
central white circles mask the pixels that are asso¬ 
ciated with a bulge and/or bar and are therefore 
not included in the fits. The velocity fields are 
given in units of kilometers per second. The ve¬ 
locity and residual velocity fields display the pix¬ 
els with measurement uncertainties less than 25 
km/s. The suppressed pixels do not strongly influ¬ 
ence the fitted parameters. The spiral effect plot 
is composed of one frame showing the response of 
the fitted i value to the rotation of the image and 
another showing how (j) changes with the rotation. 
In the asymmetry profiles, asterisks mark the A 
values for the residual velocity fields. 

Since photometric data is often more extensive 
than kinematic data, some of the images in the fol¬ 
lowing figures may appear unequal. As described 
in the Appendix, the data used by our fitting rou¬ 
tine is different for photometric and kinematic fits. 
However, we utilize the maximum amount of infor¬ 
mation in both cases. To test the impact of the dif¬ 
fering radial extents of the data, we fit the photo¬ 
metric data over the same radial range as the kine¬ 
matic data. There were no significant changes in 
the average differences between photometric and 
kinematic inclinations or position angles. For the 
remainder of the paper, the results given are de¬ 
rived using the maximum amount of information 
for the data set. 
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5.1. ESO 439G20 & ESO 322G45 

The photometric images of ESO 439G20 (Fig¬ 
ure 5a I) and ESO 322G45 (Figure 5b I) appear to 
be free of strong bars and spiral structure. Their 
respective velocity fields also appear to be those of 
simple inclined axisymmetric disks (Figures 5a II 
and 5b II). The naive assumption is that these two 
“nice” galaxies should show little asymmetry and 
should have good agreement between their photo¬ 
metric and kinematic projection angles. 

For ESO 439G20, the photometric and kine¬ 
matic i difference is 3.5° and the </> difference is 
1.6°. The lack of strong nonaxisymmetric struc¬ 
ture is evident from the photometric residual (Fig¬ 
ure 5a III) and the plot illustrating the spiral effect 
(Figure 5a V). The asymmetry profile (Figure 5a 
VI) shows little influence of bars or spirals. The 
last A value is due to the small blue and red spots 
on either side of the minor axis on the right side 
of Figure 5a IV. 

ESO 322G45 has an i difference of 4.7° and 
a (p difference of 4.5°. Figures 5b III and 5b V 
are evidence that there is little asymmetry in the 
photometric image. Again, the asymmetry profile 
(Figure 5b VI) is fairly featureless. The last two A 
values may be the effect of a weak spiral arm that 
appears in the photometric residual image « 20 
arcsec from the center. It is surprising that such a 
featureless galaxy has large i and </> differences, but 
they cannot be symptoms of intrinsic ellipticity in 
the disk, for the following reasons: An elliptical 
disk viewed in projection should generally give rise 
to the type of asymmetry profile discussed in §4.3, 
which is not seen in Fig. 5b VI. The special case 
of projection about one of the disk’s principal axes 
could mask this signature, but then we should not 
expect a significant A^, which we also find. Thus, 
we are confident that the differences in this case 
are not due to an intrinsically elliptic disk. 

These two galaxies serve to illustrate the point 
that just because a galaxy’s image or velocity field 
appears smooth or simple does not imply accurate 
estimation of projection angles. 

5.2. ESO 267G29 & ESO 569G17 

In contrast to the previous examples, we 
present two galaxies whose photometric images 
(Figures 5c I and 5d I) show strong bars and spi¬ 
ral structure. Their velocity fields (Figures 5c II 


and 5d II) contain distorted isovelocity contours, 
evidence of strong nonaxisymmetric structure. In 
these cases, it is not unreasonable to believe that 
estimation of projection angles will be problem¬ 
atic. 

ESO 267G29 has a photometric/kinematic i dif¬ 
ference of 20° and a (p difference of 15°. The strong 
open spiral can be clearly seen in Figure 5c III 
and the effect of the central bar is obvious from 
the TO = 1 residual pattern in central region of 
Figure 5c IV. The asymmetry profile (Figure 5c 
VI) clearly shows the effect of the central bar with 
the rise in A at large radius due most likely to the 
strong spiral also present. 

To again contrast, the differences between pho¬ 
tometric and kinematic i and p values for ESO 
569G17 are < 1°. A dusty multi-armed spiral pat¬ 
tern can be seen in Figure 5d III. As before, the 
effect of the central bar is evident in the residual 
velocities (Figure 5d IV). Like ESO 267G29, the 
asymmetry profile in Figure 5d VI shows the im¬ 
pact of a central bar and some effect from spiral 
arms. 

It is not surprising that if the central region of 
the velocity field of ESO 267G29 is ignored, as it 
is in the photometry, the differences between pho¬ 
tometric and kinematic i and p drop to w 10°. In 
contrast, ignoring the central region of the veloc¬ 
ity field of ESO 569G17 increases the discrepancy 
between photometric and kinematic projection an¬ 
gles by w 1°. So, even though photometric images 
and kinematic maps are complicated, the projec¬ 
tion angles derived from them can be consistent. 
These two examples also suggest that the details of 
the spiral structure present can significantly alter 
the photometric/kinematic discrepancy. 

5.3. ESO 215G39 

As our final example, we present a “typical” 
spiral galaxy. The photometric/kinematic i dif¬ 
ference is 1.9° and the p difference is 2.3°. The 
photometry and residuals (Figures 5e I and 5e II) 
show moderate spiral structure and the velocity 
field (Figure 5e II) is fairly simple. The velocity 
residual field (Figure 5e IV) shows some evidence 
of elliptical orbits in the central region. Again, 
the asymmetry profile (Figure 5e VI) contains ev¬ 
idence of either a weak central bar or spirals near 
the center. Interestingly, the magnitudes of these 
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photometric/kinematic differences are the same as 
the predicted differences due to the spiral effect for 
this galaxy. 

6. Discussion & Conclusions 

We have developed and tested a global fitting 
method with a novel minimization technique to 
obtain disk galaxy inclinations and position angles 
from either photometric or kinematic data. Since 
we fit a thin axisymmetric disk model to all the 
data at once, we determine globally optimal values 
for the geometric parameters, e and (j). This con¬ 
trasts with the IRAF ellipse utility which returns 
independent e and (j) values for each of many el¬ 
lipses, and cannot be constrained to find the best 
common values. 

We map the Xr surface to find the best-fit in¬ 
clination and position angle and to determine the 
statistical uncertainty for each fit. Due to the large 
number of data points (typically > 10^), the sta¬ 
tistical uncertainty is very small. The i and (j) 
values determined using our method differ only 
slightly from those previously estimated for the 
same 74 spiral galaxies by Palunas & Williams 
(2000). While we have held the rotation cen¬ 
ter fixed at the location of brightest continuum, 
changes in derived quantities are generally small 
when we allow the kinematic center to move from 
this point. 

A major result of this study is that systematics 
dominate photometrically determined inclination 
and position angle uncertainties. We estimate that 
systematic uncertainties due to nonaxisymmetric 
structure such as spirals and/or bars are, on av¬ 
erage, w 4° for inclination estimates and « 5° 
for position angle estimates for galaxies in this 
sample. Again, we point out that the complete 
systematic error is unknown, and these estimates 
must be taken as lower limits. It is reassuring 
that the position angle uncertainties tend to be 
smaller for more edge-on galaxies and that later 
type spirals generally have larger systematic incli¬ 
nation and position angle uncertainties. It seems 
that the only way to reduce these systematic errors 
would be to model the nonaxisymmetric structure 
in each galaxy. 

It is interesting to note that the average dif¬ 
ference in inclination and position angle derived 
from the photometric and kinematic data turn 


out to be roughly the same size as our estimates 
for uncertainties due to nonaxisymmetric struc¬ 
ture. Since these differences are also similar to 
those expected from intrinsically elliptical disks, 
they reinforce the concern raised by Franx & de 
Zeeuw (1992) and Rix & Zaritsky (1995) that spi¬ 
ral structure can be confused with intrinsic disk 
ellipticity. As discussed in §5, the (non)existence 
of spirals and/or bars does not necessarily imply 
(small) large differences between photometric and 
kinematic projection angles. We also find no clear 
evidence for the velocity field distortions expected 
for intrinsically elliptic disks. 

To conclude, we have found a difference be¬ 
tween galaxy projection angles determined from 
photometric and kinematic data for 74 galaxies. 
We have also found that nonaxisymmetric struc¬ 
tures, such as spirals and bars, impact the galax¬ 
ies’ photometric projection angles to about the 
same extent as the differences between photomet¬ 
ric and kinematic projection angles. Further, we 
see no direct evidence for intrinisically elliptical 
disks in any of our galaxies’ residual velocity fields. 
In light of these points, it seems unnecessary to ex¬ 
plain photometric/kinematic projection angle dif¬ 
ferences with intrinsic disk ellipticities. 
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A. Procedure to Estimate Projection Angles 

The principle of our method is to use as much information as possible from a photometric or kinematic 
image to estimate the two angles of projection. If a galaxy were a razor-thin, axisymmetric disk and the 
orbits were circular, then the surface brightness or in-plane orbital speed would be constant on a nested set 
of concentric, co-axial, similar ellipses. We utilize the entire velocity field, from the center of the galaxy to 
beyond its edge, but our photometric fits are to only those parts of the image not strongly contaminated by 
bulge and/or bar components. 

Our input data are a set of data values {Di}, and their associated uncertainties, {<7^}, for each of the N 
pixels in the image. For photometric data, on the one hand, the {Di\ are pixel intensity values, for which 
the principal source of noise is photon counting. We estimate the sky brightness from part of the image 
away from the galaxy and subtract its value before beginning the fit. For kinematic data, on the other 
hand, the {Di] are line-of-sight velocity values and the associated uncertainties result from the fitting of the 
Fabry-Perot data cube. 

We choose the center to have the coordinates {xc,yc) of the pixel that has the largest intensity for 
a photometric image, or the largest continuum intensity for the velocity map. We generally hold these 
coordinates fixed for the entire fitting procedure, although our method can allow the center to move if 
desired. 

We fit a non-parametric profile that is characterized by a tabulated set of values {Mk}, for a number K 
of semi-major axes. Each Mk has the same value around an ellipse, and all ellipses have a common center 
and ellipticity, e, and major-axis position angle, </>, embodying our assumption that the disk is thin and 
axisymmetric. For photometric data, the {Mk} are isophote intensities, while they are circular velocities for 
a kinematic map. Note such a tabular form makes no assmuption about the radial profile of either the light or 
the orbital speed. We derive a model prediction at a general location in the image by finding the semi-major 
axis of the ellipse that passes through the point in question, and then interpolate for the predicted pixel 
value from the {Mk}. 

For any values of the projection parameters (e, (j)), we form 
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- Y.k=l'^k,iMk 
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where the Wk,i are the weight factors that describe the interpolation scheme, and v is the number of degrees 
of freedom in the fit. Since Xr is linear in the M^s, the MkS that minimize Xr can be determined by setting 
the derivative of eq. (Al) with respect to each Mj equal to zero and solving the resulting set simultaneous 
equations. 
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The resulting Mk values are used to determine Xr ■ Oar implementation uses linear interpolation to determine 
the Wk,i factors, but more complicated schemes would be straightforward to implement. Any pixels that lie 
outside the final fitting ellipse have their model values linearly extrapolated from the values of the last two 
ellipses and their weights smoothly decreased depending on their distance from the last ellipse. Such pixels 
have little impact on the best-fit parameters. 

We use slightly different methods to determine which pixels will be used in the fitting procedure for the 
photometric and kinematic data. In both cases, we create a pixel list only once which is used throughout 
the fitting procedure. With photometric data, we exclude pixels contaminated by foreground stars or cosmic 
rays. We also exclude those pixels inside a circular aperture containing a central bulge or bar, and set the 
smallest semi-major axis of the K ellipses equal to the radius of this inner aperture. The outer limit of the 
pixel list is defined by an outer elliptical aperture with semi-major axis 10% larger than that of the outermost 
fitting ellipse, using an initial estimate of the ellipticity. The pixel list for the kinematic data contains all 
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pixels that have continuum intensity values > 0.5% the central continuum intensity. We add, in quadrature, 
a fixed velocity dispersion of 7 km s“^ to the estimated velocity uncertainty in the line fit. Some anomalous 
velocities remain in the velocity map, possibly caused by intervening gas not associated with the galaxy of 
interest or from Voigt profile fits near the velocity limits of the Fabry-Perot data cube; we eliminate such 
values by applying a simple cut to exclude pixels with velocities more than 500 km/s away from the data 
value at the fixed center. 

We vary the two projection angles, e and (j), recomputing the at each iteration to determine their 

values that minimize Xr ■ We use a downhill simplex method, specifically the amoeba routine described in 
Press et al. (1992), which requires only function evaluations (no derivatives) and initial estimates for e 
and 4>. Following the recommendation of Press et al. (1992), the e and </> estimates returned are re-used as 
initial guesses for up to 5 subsequent calls to amoeba to converge to the Xr minimum. We have verified the 
minimum location by inspecting the Xr surface surrounding the best-fit values. Projections of the contour 
where Xr = Xrmin + 1-0/j^ onto the e and <j) axes provide 1-a statistical uncertainty measurements. 

Figure 6 shows the results of photometric fitting for one galaxy in our sample (ESO 215G39). The model 
isointensity values are shown as filled circles superimposed on data values. The solid lines encompass 68% of 
the data values. The statistical uncertainties of the parameters (dot values) are smaller than the symbols. 
The rotation curve with 1-a parameter uncertainties resulting from the kinematic fit for the same galaxy is 
shown in Figure 7. 
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Fig. 1.— a) Inclinations determined from pho- 
tometric data (iphot) compared with those deter¬ 
mined from kinematic data (tvei)- b) Position 
angles determined from photometric data ((/iphot) 
compared to those determined from kinematic 
data {4>ve\)- The extreme outlier is discussed in 
the text. The solid lines in both frames are lines 
of equality. 


< 



Tphoi 


Fig. 2.— a) The difference between inclination an¬ 
gles determined from photometric data and kine¬ 
matic data plotted against the photometric in¬ 
clination. b) The difference between position 
angles determined photometrically and kinemat¬ 
ically versus the photometric inclination. Note 
that the scatter shrinks as the inclination in¬ 
creases. The extreme outlier is discussed in the 
text. 
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Fig. 3.— a) Histogram of the distance between 
the photometric center and the fitted kinematic 
center. The difference is < 2."7 for the majority 
(90%) of the galaxies, b) Histogram of the change 
in estimated inclination between the fixed-center 
and floating-center fits. 95% of the galaxies have 
Ai < 7°. c) Histogram of the change in position 
angle between the fixed-center and floating-center 
fits. None of the galaxies has Acp > 4°. 
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Fig. 4.— Systematic uncertainties in (a) photo¬ 
metric inclination, (b) photometric position an¬ 
gle. The numbers denote the T-type of the galaxy. 
Note that the scatter decreases with increasing 

^phot • 


Fig. 5.— (I) /-band image of the galaxy. The 
white circles show the locations of removed fore¬ 
ground stars. The units are mag/arcsec^. (H) 
The line-of-sight velocity field from Fabry-Perot 
interferometry. The values listed are heliocentric 
velocities in km/s. (HI) Photometric residuals us¬ 
ing the best-fit Zphot and ^phot- The units are 
mag/arcsec^. The ellipse marks the location of 
the last ellipse used in the fitting routine. The 
central white circle masks the pixels that are asso¬ 
ciated with a bulge and/or bar. (IV) The residual 
velocity field using the best-fit Zvei and 4>vei- As 
explained in the text (see §5), only those pixels 
with measurement uncertainties < 25 km/s are 
shown. (V) The variation of the fitted i and (p 
as the image is rotated as described in §4.2 (VI) 
The radial variation of A to reveal the asymmetry 
profile. Discussion for each galaxy is given in the 
text. 

Fig. 5.— Continued 
Fig. 5.— Continued 
Fig. 5.— Continued 
Fig. 5.— Continued 
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Fig. 6.— Model isophote values {Mk} (filled 
circles) versus semi-major axis for the galaxy 
ESO 215g39. Some indicative (but not all) sky- 
subtracted data values are plotted. The dashed 
lines contain 68% of the data values and the 1- 
a uncertainties for the model values are much 
smaller than the filled circles. 
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Fig. 7.— Estimated rotation curve for the galaxy 
ESO 215g39. The crosses with error bars are the 
{Mfc} values returned from the kinematic fitting 
routine together with their 1-a statistical uncer¬ 
tainties. 
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